Low temperature decoherence by electron-electron interactions: 
Role of quantum fluctuations 
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We derive a general expression for the conductivity of a disordered conductor with electron- 
electron interactions (treated within the standard model) and evaluate the weak localization correc- 
tion (5awi employing no approximations beyond the accuracy of the definition of Sa^i. Our analysis 
applies to all orders in the interaction and extends our previous calculation by explicitly taking into 
account quantum fluctuations around the classical paths for interacting electrons (pre-exponent). 
We specifically address the most interesting low temperature limit and demonstrate that such fiuc- 
tuations can only be important in the perturbative regime of short times while they are practically 
irrelevant for the Cooperon dynamics at longer times. We fully confirm our conclusion about the 
existence of interaction-induced decoherence of electrons at zero temperature for the problem in 
question. We also demonstrate irrelevance of a perturbative calculation by Aleiner et al. (AAV) [J. 
Low Temp. Phys. 126, 1377 (2002)] and refute AAV's critique of our earlier analysis. 



I. INTRODUCTION 



The electron decoherence time in disordered conduc- 
tors saturates to a finite value at low temperatures 0. 
In Ref. we offered an explanation of this effect at- 
tributing it to electron-electron interactions. This expla- 
nation was supported by a detailed nonperturbative cal- 
culation and the results 0^ are in a good agreement 
with experimental findings 1]. Subsequently Aleiner et 
al. (AAG) developed an alternative ~ perturbative in the 
interaction - calculation Q and claimed that the results 
of the latter {i) contradict to our results 0, 0] and (m) 
yield zero electron dephasing rate at T ^ 0. AAG also 
claimed (Sec. 6.1 and 6.2 of H) that they have found a 
"mistake" in our calculation Q. 

All these claims have been carefully analyzed and 
demonstrated to be in error j^. We have argued that 
(a) the perturbative calculation Q yields ambiguous re- 
sults and, hence, is useless for the problem of electron 
dephasing at low temperatures, (b) even if one adopts 
the perturbative strategy Q one recovers a finite elec- 
tron dephasing time at T ^ 0, (c) on a perturbative 
level our results Q do agree with those of Ref. and 
(d) the claim ^ about "missing diagrams" in our calcu- 
lation is wrong. This discussion and further arguments 
on low temperature dephasing due to electron-electron 
interactions were reviewed in Ref. 0. 

Recently Aleiner et al. (AAV) 7] have made another 
attempt to challenge our results and conclusions. One 
of the goals of the present paper is to analyze and re- 
fute these new criticisms. In particular, we will point out 
that the replacement of "the density matrix by its Wigner 
transform", eq. (5) of ^7], claimed by AAV as "the main 
source of the mistake" was not performed at all in our 
derivation of the effective action, cf. eq. (43) of Ref. 3J. 
We will also analyze and prove irrelevant another sugges- 
tion of AAV, that the term 5r in the effective action Q 
should contain an imaginary part which - according to 



AAV - would be responsible for the well-known cancella- 
tion of "coth" and "tanh" in the first order perturbation 
theory at T ^ 0. 

Since AAV's critique of our calculation is essentially re- 
stricted to these two claims |£| , the above observations 
are already sufficient to conclude this discussion. How- 
ever, taking into account fundamental importance of the 
issue, we have performed an additional analysis aimed to 
construct a complete solution of the problem. Our main 
goal is to evaluate the weak localization (WL) correction 
to the conductivity in the presence of interactions making 
no approximations beyond the accuracy of the definition 
of this quantity. This solution is worked out below and 
our main result is presented in eqs. (|40|l - (|42[l . 

Our paper is organized as follows. After brief remarks 
on physics and experiment (Section 2) we demonstrate 
(in Section 3) that AAV's perturbative calculation Q is 
unsuitable for the problem of quantum decoherence of 
electrons at low T. In Section 4 we briefly recollect our 
earlier results and, making use of general arguments, ex- 
plain why our path integral analysis 0, 0, U is sufficient 
for the problem in question. The main results of this pa- 
per are presented in Section 5. We first derive a formally 
exact expression for the conductivity of an arbitrary dis- 
ordered conductor in the presence of electron-electron in- 
teractions (Sec. 5A) and then use it to explicitly eval- 
uate the WL correction to all orders in the interaction 
(Sec. 5B and 5C). This analysis confirms our previous 
results 0, Q and extends them by fully accounting for 
quantum fluctuations around the classical paths for in- 
teracting electrons. We compare our results with those of 
other authors and present further discussion in Sec. 5D 
and 5E. In Section 6 we specifically address and refute 
AAV's critique T] of our earlier calculation. A brief sum- 
mary is presented in Section 7. Some technical details are 
relegated to Appendix. 
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II. PHYSICS AND EXPERIMENT 

In Ref. laAAVpointed out that our results and con- 
clusions UlS 1^1 M "physically inconsistent" with 
the qualitative arguments [J against quantum dephas- 
ing in the zero temperature limit. According to these 
authors quantum decoherence would be impossible in a 
subsystem of any interacting quantum system provided 
the latter is close to equilibrium at T — > 0. At this point 
we note that the arguments 4] contradict not only to our 
conclusions but also to numerous results for various other 
models - including exactly solvable ones - where quan- 
tum decoherence of one degree of freedom is obtained as 
a result of its interaction with others even at T ^ 0, see, 
e.g., Hence, no general proof can be constructed 

which would rule out interaction-induced quantum deco- 
herence at T ^ 0. The conclusion about the presence 
or absence of zero temperature dephasing in any given 
model can only be obtained from a detailed calculation, 
not on the basis of "general arguments" . 

AAV also claimed ItI that there exists an "overwhelm- 
ing experimental evidence" against our statements. By 
means of a detailed comparison with the experimental 
data we have demonstrated that our predictions arc in a 
good quantitative ^reement with numerous experiments 
(see, e.g., Refs. HI SB ^2) including those (see Section 
5 of Ref. Bl) where for some samples no saturation of 
Tip was claimed down to T ^ 50 mK. Therefore, a bulk 
piece of existing experimental results - including very re- 
cent ones ^3, - clearly supports the conclusion about 
the presence of interaction-induced dephasing at low T 
rather than argues against it. 



III. INSUFFICIENCY OF AAV'S 
PERTURBATIVE APPROACH 

In Refs. 0, B| we have already explained in details 
why perturbative in the interaction techniques - espe- 
cially if combined with the Golden rule approximation 
- are insufficient for the problem in question. However, 
since AAV still keep using such techniques, we will briefly 
repeat our arguments adopting them to the calculation 



repe 
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In Ref. 7] AAV have related the dephasing time for 
an electron in a disordered conductor to the self-energy 
for the Cooperon. Expressing the Cooperon in the form 



T ^ 0, since according to Eq. (1) of Ref. in this limit 
one should set e — > 0. Another drawback of the definition 
(0) is that the self-energy is evaluated only in the limit 
(5 = 0. However, the dependence of S on the momentum 
Q in eq. |^ cannot be a priori neglected. 

In addition to the above approximations, AAV sug- 
gested to replace the exact self-energy in Eq. ^ by the 
result of the first order perturbation theory in the inter- 
action S'-'^' (cj, Q, e) : 



rAAV 



= s(i)(w = o,g = 0,e = 0). 



(3) 



This suggestion constitutes their major mistake. Quite 
obviously, it is not possible to recover the unknown func- 
tion (in our case E) if one only evaluates the first order 
term (S^^^) of its Taylor expansion. 

In order to illustrate this point it suffices to consider 
the following simple example. Assume, for instance, that 
the "Cooperon" C{t) depends weakly on the coordinates 
and neglect this dependence. We also assume that in the 
presence of interactions this "Cooperon" decays in time 
as 



C(t) = 6i(i)(l + at)e-("+'5'^)*. 



(4) 



where a and /? are proportional to the interaction 
strength, i.e Q reduces to 6{t) in the absence of interac- 
tions. After the Fourier transformation of Q one readily 
finds 



C{lu) 
S(c.) 



1 



—iuj + Tt{uj) 

(a + PTf - iLuPT 



2a + f3T- ioj 

Combining with the definition Q we obtain 

1 _ a + (3T 
~ ~ 1 



(5) 



(6) 



'"P a+(3T 

i.e. at T ^ one arrives at a non-zero dephasing rate 

Let us now evaluate for the same example following 
the approach of AAV. For this purpose we expand the 
exact expression for Ti^uj) (jSJ in powers of interaction 
and get 



C(w, (5,e) 



(1) 



they have defined the dephasing time at T = as follows 



= I](c^ = 0,Q = 0,e), 



(2) 



see the unnumbered equation after eq. (2) of This 
definition is ambiguous since the self-energy is a function 
of e whereas in our problem is a function of temper- 
ature T but not of e. The ambiguity disappears only at 



t(Lo) =I3T + 



(8) 



Keeping only the first order contribution to the self- 
energy, from AAV's eq. (jSJ we find 



1 



rAAV 



f3T. 



(9) 



This result differs drastically from the exact one © at 
sufficiently low temperatures. In particular, at T — > the 
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dephasing rate ® vanishes, while the exact expression 
approaches a non-zero (linear in the interaction) value 
(|7| . Furthermore, all higher order terms in the expansion 
© do not vanish at T and, moreover, diverge at 
small oj. 

The reason for the failure of AAV's perturbative ap- 
proach is, of course, obvious from eq. lO: An expansion 
of this expression in powers of a and (3 is only justified 
for uj ^ a + PT, i.e. in the limit of high frequencies or 
short times. We, in contrast, are interested in the oppo- 
site limit of small frequencies or long times. Thus, AAV 
simply missed the low temperature contribution to the 
dephasing rate by inadequately extending their perturba- 
tive expansion of Tt{uj,Q,e) to low frequencies, whereas 
it can be applied at high frequencies only. 

Our simple example, eq. (0), also illustrates irrelevance 
of the claim that our approach "is equivalent to calcu- 
lating only a single contribution S*^'') to the self-energy 
and using the conventional Dyson equation". AAV ar- 
rived at this conclusion simply by observing the combi- 
nation S^*^^ — E^'^"^) ~ PT in their first order perturba- 
tive result (cf. eqs. (3,4) in 0) and comparing it with 
our result, which contains not only /3T but also the T- 
independent contribution. The only - superficial - rea- 
son for Aleiner et al. to perform this comparison and to 
qualify our results as "purely perturbative" is that our 
expression for the dephasing rate "is proportional to the 
first power of the fluctuation propagator" . We hope it 
should be sufficiently clear from eq. © that linear in the 
interaction expression jT)) is non-perturbative and, hence, 
it cannot be obtained by a simple expansion JH)) of the 
self-energy S(w) in the interaction. 

In order to avoid misunderstandings let us empha- 
size that the above example (0J is not meant to be an 
explicit solution for the problem with electron-electron 
interactions. This solution will be worked out below in 
Sec. 5. Eq. ^ is just an illustration of one of the draw- 
backs of the perturbative approach to the problem of 
quantum dephasing. In the problem with disorder and 
electron-electron interactions the situation turns out to 
be by far more complicated. For instance, already the 
first order result diverges both with time and at large 
frequencies in Id and 2d systems, see, e.g., eq. (70) of 
Ref. 1^]. 



X e'i 



^{iSolpi,xi]-iSolP2,X2]-iSRlpi,xi,p2,X2]-Silxi 



^2])(10) 



Here So[p, x] = dt' [px — / 2m — Uirn-p{x)\ is the action 
of a nonintcracting electron in a disordered potential of 
impurities and the terms iS^ and Si describe the effect 
of the bath (formed by all the electrons) on the motion 
of a single electron. This form of the effective action 
for an interacting particle is standard in the Feynman- 
Vernon theory of influence functionals The path 

integral (|10|l is evaluated within a semiclassical approx- 
imation controlled by the parameter kpl 3> 1. For a 
d-dimensional conductor we find 



c{t,x = 0) = e{t) 



Ad{t) 



-hit) 



(11) 



where D is the diffusion coefhcient, Ad{t) is pre-exponent 
{Ad{t) = I without interactions), and the function fd(t) 
is obtained by evaluating the influence functional on pairs 
of classical time-reversed paths. The result can be writ- 
ten in the form Hi 



fdit)=at + 6.fdiT,t), 



(12) 



where the function 5fd{T,t) is oc 

Tt3/2 for d = f and 
cx TtlnTt for d = 2 in the Hmit > 1, while it is oc 
\/t\nt for d = 1 and cx Int in the opposite limit Tt -^1. 
The linear dependence 5fd{T,t) cx Tt strictly applies for 
d = 3. 

It is also instructive to explicitly indicate the depen- 
dence of our results on the Plank's constant fi. At T ^ 
we have 



HCd{t)/cf\t)] 



S'(=i)(i) 



HAd{n,t% (13) 



where S'^'^^^ ~ at is the classical (?i-independent) action 
on time-reversed saddle point paths (a is ?i-independent 
and a = a/K) and TiXnAd represents the quantum cor- 
rection to the classical action. This quantum correction 
can only be important if S'^'^'-' is small as compared to 
i.e. at times t <^ Ti/a. The perturbative approach 
Q applies only in this limit. On the other hand, for 
nonzero t the quantum correction TihiAd should vanish 
at ?i ^ 0. Already because of this reason it cannot cancel 
the classical part of the action S'^'^'^ > h. It is therefore 
sufficient to evaluate the ?i-independent part of — z5r — Si 
on pairs of classical time-reversed paths and obtain the 
IV. PATH INTEGRAL ANALYSIS: EXPONENT dephasing time from the condition ^(=i)(t^) ~ ?i. Fur- 
thermore, it turns out that S'r vanishes for such paths 
implying that S'r can only contribute to the pre-exponent 
but not to S'f'^''. Hence, S'r, is irrelevant for and the 
idea that quantum fluctuations around the time-reversed 
classical paths generate the S'^-dependent contribution 
to the classical action can be rejected on general grounds 
without any calculation p^ . 

This is the logics behind our saddle point analysis 0,0 
which AAV attempted to challenge. The only way to 
support the AAV's arguments is to prove that at T 
the term S'r provides a contribution to Ad proportional 



An important advantage of our path integral approach 
is the possibility to describe the long-time behavior of the 
Cooperon with exponential accuracy, which is sufficient 
for the problem in question. This approach is free from 
ambiguities inherent to the perturbation theory ^ 2\ ■ 

We define the kernel of the evolution on the Keldysh 
contour in terms of the path integral |H 

/ Vpx I Vxi I Vp2 I VX2 
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to exp{at/h) which grows exponentially with time and 
diverges in the classical limit ?i — > [in Ref. j3l such 
a contribution was claimed to be provided by the terms 
E^'^"^)]. Below we will explicitly evaluate not only 5^^^'^ 
but also the quantum correction to the effective action in 
all orders in the interaction. We will demonstrate that 
- in contrast to the AAV's claims - the 5R-dependent 
pre-exponent cannot grow at sufficiently long times and, 
hence, in no way can compensate an exponentially de- 
caying contribution from the 5'i-terms. 



V. DECOHERENCE BY INTERACTIONS 

A. Exact results 

In the beginning we will closely follow the analysis of 
Ref. (3 where the interested reader can find further de- 
tails. We start from the general quantum mechanical 
expression for the linear conductivity a which can be 
written in the form 



dt' (tr (ji^)Uiit, t')[x, pv{t')]U2{t\ t) 



3ih 

Here the current density operator is defined via 
he 



(14) 



{Xi\j{x)\x2) 



[^a:i ,xSx2,x ^xi ,x^ X2 .x 

] (15) 



and implies averaging over the fiuctuating quantum 
fields V'^ and which mediate Coulomb interaction 
between electrons. In eq. (|14|) and below we implicitly 
assume averaging over the coordinate x where the current 
density is evaluated. 

The evolution operators Ui,2 in (|14(l are 



lJi,2{t,t') = Texp 



dTH,^2{T) 



(16) 



where 



Hi{t) = i?o - /i - eV+it) - -[1 - 2pvit)]eV-{t), 
H2{t) ^Ho-P- eV+{t) + ^eV-{t)[l - 2pv{t)], (17) 

and Hq = /2m + C/imp(i) is the Hamiltonian of a non- 
interacting electron. The density matrix pv{t') obeys the 
non-linear equation 

^ i^o - eV^ , Pv] - ^{1 - Pv)eV^ Pv 



^PveV (l-pv)- 



(18) 



Averaging of pv{t) over yields the exact single 
electron density matrix in the presence of disorder and 
Coulomb interactions. The next steps in Refs. 0, 



were to express the kernels of the operators C/1.2 in terms 
of the path integrals and average their product over the 
fiuctuating fields V^. After that one arrives at the path 
integral ((Tn|l . 

At this point we depart from the analysis of Refs. [1,1^ . 
We will postpone using the path integrals and continue 
exact manipulations with the operators. We first note 
that the solution of eq. H18|) with the initial condition 

py{0) = po = [1 + e(^°"^'/^]"i can be expressed in the 
following exact form 



pvit) = [l + U2(t,0)e(^''-'')/^ui(0,t) 



where we have defined 

ui,2it,t') = Texp 

and 



Hi.2 = Hoit') -p- eV+{t') T eV'{t')/2. 



(19) 



(20) 



(21) 



We then observe that the operators U1.2 satisfy the 
Schrodinger equation 

ih^UiitX) - Hiit)Uiit,t'), Ui{t',t') = i, 

ot 

in^U2{t',t) = U2{t\t)H2{t), U2{t',t') = l. (22) 

The solutions of (|22|l can be found exactly. They are 

Ui{t,t') = [l-pvit)]uiit,t')[l-pvit')]-\ 

U2{t\ t)^[l- pv{t')Vu2{t',t)[l - pv{t)]. (23) 

Combining these expressions with eq. H19|l one can 
rewrite the operators H23(l in the following identical form 

J7i(i,0 = [{\-pv{t'))u^{t',t) + pv{t')u2{t',t)]-\ 
U2{t',t) - [U2{t,t'){l - pv{t'))+uiit,t')pv{t')r\{2A) 

Eqs. (|24|) is our key technical result. Together with eqs. 
(|14I15|I this result provides a formally exact expression 
for the linear conductivity of an arbitrary disordered con- 
ductor in the presence of electron-electron interactions. 
All the diagrams of the perturbation theory in all orders 
in the interaction are fully contained in the above ex- 
pressions and can be recovered by expanding (|24|l in 
with subsequent averaging over these fields. For instance, 
expanding H24I) to the second order in V^, after aver- 
aging one arrives at contributions of the type {V~^V~^) 
and {V^{l — 2p)V^) which yield respectively "coth" and 
"tanh" terms in the perturbation theory |j, 01 ■ 

Thus, making no approximations we have demon- 
strated that the time evolution of the single electron den- 
sity matrix in the presence of interactions is determined 
by the operators ui^2{t,t') (|2U|I which do not contain the 
density matrix pv at all. We also note that the operators 
Ui,2 depend on the density matrix pv taken at the 
initial time t' only but not at later time moments. Below 
we will make use of these features and evaluate the WL 
correction to conductivity to all orders in the interaction. 
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B. Quasiclassics 

Let us first prepare tlie main building blocks of our 
calculation of the WL correction. This calculation will 
then be completed in Sec. 5C. 

To begin with, we notice that in eqs. 124() relative 
contributions of the evolution operators ^1,2(^1^') are 
"weighted" by the factors l—pv{t') and pv{t')- In the low 
temperature limit almost at any electron energy one of 
these factors dominates over the other p^. Hence, one 
of the two operators ui or U2 in H24|l can be neglected 
except if the eigenvalues of pv{t') are close to 1/2. Con- 
sider, e.g., small eigenvalues of this operator. This situ- 
ation describes electrons with energies above the Fermi 
level. In this case terms containing pv{t') hi (|24|l can be 
neglected and we arrive at the following contribution to 
{xi\Ui{t, t')[x,pvit')]U2{t',t)\x2): 

J J dzidz2{xi\ui{t, t')\zi)(z2\u2{t' ,t)\x2) 

x{zi\[x,pv{t')]\z2). (25) 

Let us express the matrix elements of the operators ui,2 
via the path integrals 

{yi,2\ui.2{t,t')M = y I?x(T)e^^i-^[^Ml , (26) 
1/1,2 

where 5i.2 are the exact actions pertaining to the Hamil- 
tonians H21|) 

Si^2^So + eJ^^ dT[V+iT,x{T))±V~iT,xiT))/2], (27) 

So - £ dr (^a^^ _ C/i^p(x(r))^ . (28) 

We emphasize again that eqs. (|26l28f) are exact and they 
do not contain the electron density matrix pv at all. We 
should now (a) evaluate the path integrals H2()|l and then 
(b) average the combination (|25() over the fluctuating 
fields V^. 

Evaluation of the matrix elements. Let us make use of 
the fact that the WL correction to conductivity is defined 
within the accuracy kpl ^ 1. This inequality is usually 
well satisfied in disordered metallic conductors. Hence, 
we can evaluate the matrix elements (|26|l quasiclassically. 
Since the actions l|27|) do not depend on pv one can con- 
veniently employ a regular expansion of 6*1 _2 in powers of 
h. The path integrals are then easily evaluated and we 
arrive at the well-known van Vleck formula 

{Xl,2\ui,2{t,t')\zi^2) 

= y V A det- ^ e« ^i.^ , (29) 

^ V \27r?i J dxi 2dzi 2 ^ ' 

where s'(l^.^ = (t, i', a;i(2) , 21(2)) = S'i(2) [5i„(2„)] and 



^1,2 are the exact least action paths obeying the equa- 
tions 

^Si(2)[xin(2n),y^]/ Sxin(2n) =0 (30) 

with the boundary conditions iin(2n)(i') = ^1(2.) and 
Xin(2n){i) = -^1(2) • In general there exist several or even 
many different classical paths satisfying the above condi- 
tions. Here and below the index n labels all such paths. 

Clearly, eq. (|29|l accounts not only for the saddle point 
trajectories (exponent) but also for quantum fluctuations 
around xin,2n(T) (pre-exponent) for arbitrary . It is 
also completely obvious that there is no way how the pre- 
exponent can cancel the exponent for any configuration 
of the fields . Hence, such cancellation is impossible 
also after averaging over these fields no matter what the 
details of this averaging procedure are. 

Averaging over the fluctuating fields. In general av- 
eraging of the combination 125|) over involves path 
integrals ljAl|l over these fields at all times from zero to 
t. This is because pv{t') is nonlocal in time: Accord- 
ing to eq. (|19ll it depends on times between zero and 
t'. However, with the accuracy kpl ^ 1 one can perform 
averaging in l|25|) at times smaller and larger than t' sepa- 
rately. This splitting is achieved by expressing pv{t') H19() 
via the path integrals, making use of eqs. 12t)l28|l and av- 
eraging the whole combination H25|l over at all times 
between and t. One arrives at the effective actions con- 
taining nonlocal in time contributions >S'Rj(ii, ^2) which 
vanish for all the trajectories relevant in the quasiclassi- 
cal limit kpl ^ 1 provided ti > t' , t2 < t' and vice versa. 
As a result we obtain 

dzidz2Ji2{t,t';xi^2,zi^2)izi - 22)^(21,22), (31) 
J12 = {{xi\ui{t,t')\zi){z2\u2{t\t)\x2))v. (32) 

Here we used (zi|[i;,/5]|z2) = {zi — 22)^(21,22), where 
p = {pv{t'))v is the exact equilibrium electron density 
matrix in the presence of interactions. We can also add 
that, as it was explained in Sec. 4 of Ref. J5], with the 
same accuracy kpl ^ I one can replace 

/5y(i')— Po = [l + e(^°-'^'/^]-^ (33) 

already before averaging over . After that the factor- 
ization (|31|l (with p po) is, of course, an exact proce- 
dure. 

What remains is to average the product of the two 
matrix elements in (|32|l . This averaging is carried out 
in a standard manner. If the fields vary at scales 
exceeding the elastic mean free path, one can neglect 
the dependence of both the classical paths ii_2 and the 
pre-exponent in eq. H29|l on the fields V^. This approx- 
imation is sufficient for evaluation of the WL correction 
to the conductivity. 

Then averaging with the action ljA2p can be performed 
exactly. We first integrate over V'^. As both actions 
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[2|) and H27|) are linear in V"*", this integration yields 
the (5— function 



S{V {t, x) -Vo{T,X,Xin{s),X2mis)), 



where 



Vo{T,X,Xin{s),X2m{s)) ^ ~e / ds[R{s - T, Xinis) - x) 



-R{s - r, X2m(s) - 



(34) 



and the function R{t,x) is defined in (jA4|) . Due to this 
(5— function the subsequent integration over V~ also be- 
comes trivial and we obtain J12 — J2n m "^12"' where 



1 



'det- — 2_det- ° 



{2Tr?i)^ y dxidzi "~ dx2dz2 

y^Q^{iSa[xin]—iSo[i2m]-iSn[iln,X2m] — Sl[xi„,X2m]) (35) 

In eq. (|35|) the term Si is identical to one derived in 
Ref. Q (see eq. (55) of that paper) while the action 
is obtained from eq. (54) of Q by formally setting 
the function n(p, r) equal to zero in that formula. The 
action Six. is purely real for any pair of paths xin and X2m- 
Thus, together with the terms iSq it can only provide 
oscillations of the kernel J and in no way can compensate 
its decayJ^j™ ^ exp(— S"!/?!). As it was already discussed 
in Ref. Q and elsewhere, the action Si is real and positive 
for any pair of trajectories (except for identical ones in 
which case Si — 0). The length of electron trajectories 
in a metal always grows with time since electrons move 
with a constant velocity ~ vp- Hence, for any pair of 
time-reversed paths Xin{s) = X2m{t + t' — s) the action 
Si grows with time as well. This in turn implies that for 
such paths the kernel decays with time and vanish 
in the long time limit t — t' — > oo at any temperature 
including, of course, T = 0. 

In the above analysis we neglected the evolution op- 
erator U2 [ui) in the exact expression for C/i (i/2) (|24|l . 
This is correct at low T and for the electron energies 
above the Fermi level. Below the Fermi energy, on the 
contrary, one can drop terms containing 1 — pvit') be- 
cause in this case the eigenvalues of pv{i') are close to 
one. Then the whole analysis is repeated, one should 
only interchange the operators ui and 112. In this way 
one again arrives at eq. (|35|l (with S'r —S-n) which 
again decays as eyi'p{—Si/K). The remaining options are 
to neglect either iii or U2 in both expressions for Ui^2 
(|24|l . One again finds the contributions cx exp{~Si/fi). 
Since in all these cases Si remains the same, one con- 
cludes that if the operators ui and U2 yield comparable 
contributions (in which case the exact form of H24|l should 
be used) they will also decay as exp(— 5'i/7i) on any pair 
of time-reversed paths. Below we present an explicit cal- 
culation of the WL correction which fully confirms this 
conclusion. 



C. Pre-exponent and weak localization correction 

We are now prepared to evaluate the conductivity. As 
before, we assume that kpl ^ 1 and that the fields 
vary in space at scales exceeding the elastic mean free 
path I. In this case quasiclassical electron trajectories 
are not disturbed by interactions, and the contributions 
of the fluctuating fields add up independently. There- 
fore, we can approximately split the operators 

ui,2(t, t') ~ Mt, t')rs{t, t', V+)sit, t', ±V-/2), (36) 

where uo{t, t') is the evolution operator pertaining to the 
non-interacting Hamiltonian and 



s{t,t',V) = uoit',t)Texp 



(37) 



Within the same accuracy we can replace 



= Texp (^-- dT[Ho - eV+iT)]j . (38) 

Combining ||2ni)-(ISHI> with we obtain 

Ui{t,t') ~ uit,t\V+){{l-pvit'))s-\t,t\V-/2) 
+pv{t')s-\t,t\-V-/2)}-\ 

U2it',t) ^ {s{t,t',-V- /2)il^ pvit')) 

+sit, t\ V- /2)pv{t')Y^u{t\ t, V+). (39) 

Substituting (|39|l into eq. I|I4|I . evaluating the matrix 
elements of the operator u{t, t' , V~^) by means of the van 
Vleck formula H29|l and integrating over the fluctuating 
fields exactly as in Sec. 5B we find 



-y. 

3™ ^-^ 



dt' j dyidy2 j dzidz2 



X (Va;i - Va;2 ) Ui =2^2 2:1, X2; 2/1,^2) 

X A'l'^(t,t' ,yi,zi,xi,X2){zi - 2:2)^0(2:1,^2) 
xAr{t,t',z2,y2,xi,X2), (40) 



where 



j'wm 



(27ra)3 



'det 



925', 



in) 



dxidy] 



■det 



(m) 



X exp<j -Sa[xin] - j^Sa[x2, 



dx2dy2 
, 1 



Sl[xin,X2m] K41) 



and 



Ar^{yi\{{i-Po)s-\t,t',v-/2) 

+ PoS^^it, t', -V^/2)} ^\zi) \v-=Vo{r,x,xi„,X2,r,) > 

Ar = {z2\{rsit,t',-V-/2)il-po) 

+ sit, t', V-/2)poy'\y2)\v-=v„ir,x.M.Mj'^^) 
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As before, the paths xin and X2m satisfy the Newton 
equation (|3UII (with — 0) and the boundary condi- 
tions ii„(2m)(i') = 2/1(2) and ii„(2m)(t) = a;i(2). 

Eqs. I|4U|) - (|42|) represent the central result of this pa- 
per. They determine the linear conductivity of an arbi- 
trary disordered conductor to all orders in the electron- 
electron interaction. The above equations are based on 
the exact results H24|l and are valid in the quasiclassical 
limit kpl S> 1. We would like to emphasize that no qua- 
siclassical approximation for the electron density matrix 
was employed during our derivation and no averaging 
over impurities was performed at all. 

Let us briefly analyze eqs. 1001- EU- According to 
the standard arguments two types of classical paths a;i„ 
and i2mj identical and time-reversed ones, play an im- 
portant role in the quasiclassical limit kpl ^ 1- For a 
pair of identical paths xin{s) = X2n{s) the two actions 
in the exponent (|41|l cancel each other, the term 5*1 
vanishes identically and the matrix elements Ai^2 reduce 
to (5-functions A™2 = S{yi^2 — ^1,2) because Vb = in 
this case. In this way we recover the well known prop- 
erty that the diffuson does not decay in time even in the 
presence of interactions. 

Here we are interested in the quantum correction 
to conductivity arising from the time-reversed paths 
iin(s) — X2m{t + t' — s) . FoT any pair of such paths 
the actions Sq cancel again but the interaction term Si is 
now positive, it grows with time and yields (exponential) 
decay of the quantity J"™ H41|l in the long time limit. 
The matrix elements ^4"™ also depend on the interaction 
and on time in this case. It is obvious, however, that A"™ 
cannot grow at long times because the function Vq (|34|I is 
purely real and, hence, s{t, t' , ±Vo/2) are the unitary op- 
erators. The matrix elements of such operators can only 
oscillate provided the function Vq changes in time. Hence, 
no compensation of decaying J"™ cx exp(— S'i(t — t')/h) 
can be expected for sufficiently large t — t', and the whole 
expression under the integral over t' in eq. (|4U|) decays 
exponentially together with J"™ for any pair of time- 
reversed paths. Obviously, the matrix elements (|42|) also 
cannot grow if one formally takes the limit ?i — > 0, while 
jnm vanishes in this limit. All that implies that - in 
full agreement with our general arguments (Sec. 4) - one 
can indeed obtain the dephasing time from the condi- 
tion Si{t^) h which does not depend on the density 
matrix pv and identically reproduces our earlier results 
[2, 3, 5, 6]. 

Although in principle one can proceed further and un- 
der certain approximations evaluate the matrix elements 
for pairs of time-reversed paths, we will not do it 
here. The reason for that is obvious from the above dis- 
cussion: Particular values of are irrelevant for de- 
phasing. It suffices to observe that these matrix elements 
do not grow at long times. 



D. Relation to other results 

It is useful to compare eqs. (|40|I - H42|) with the results of 
some earlier calculations of the WL correction. Altshuler, 
Aronov and Khmelnitskii (AAK) considered the elec- 
tron dephasing by a fluctuating external field. Their re- 
sults are easily recovered from our calculation if one sets 
V~ = 0. Then one again arrives at eqs. (|40|I - H42|I with 
= (S(yi,2 — 21,2)- AAK furthermore applied their re- 
sults to the problem of interacting electrons identifying 
an external field with one produced by fluctuating elec- 
trons (the field V'^ in our analysis). In order to account 
for Pauli blocking AAK suggested a phenomenological 
procedure which amounts to keeping only the classical 
part of this field and to cutting out its quantum modes 
with frequencies uj > T (i.e. all modes at T ^ 0). This 
last step has no analogy in our calculation. 

An attempt to justify the procedure (2i| was recently 
undertaken by AAG (Jj within the framework of the first 
order perturbation theory in the interaction. Our general 
expressions, if expanded to the first order, yield eq. (|A6p 
which can be written in the form 

a = a(°)+5af^-t-,s4^^ (43) 

Here tr*-'^^ is the non-interacting contribution defined by 
the first term in (|A6I) . Jcr}^' corresponds to the terms in 
(jA6|l which contain the product V^V^ (S'l-terms), while 
5a'^^ is given by the terms containing {\ — 2po)V^ and 
(1 — 2pQ)y~V'^ (^R-terms). As we have already shown 
in Ref. |5j|, eq. ljA6|l is exactly equivalent to one derived 
by AAG In particular, it contains the combination 
coth 1^ -I- tanh ^-^^ leading to partial cancellation of the 

terms 5a^^ and <5(T^^^ 

The perturbative result l|A6p is reproduced at every 
stage of our analysis. 

• In order to obtain (|A6|I one can just evaluate the 
operators Ui^2 perturbatively starting directly from 
their definition 1161) . Substituting the result in eq. 
H14|l and replacing pv{t') po in l|14|l one arrives 
at l|X6)l . 

• Alternatively, one can also expand the exact ex- 
pressions for J7i,2 H24I) in V^. One recovers the 
same resuh 

• One can also expand approximate expressions for 
Ui^2 II39|I with the same result. 

• Finally, one can perform a perturbative expansion 
of the quasiclassical result H40|l . One should ex- 
pand J"™ to the first order in 5*1 and A"™ to 
the first order in Vq- One obtains terms propor- 
tional to the functions (V^V^) I (jASp and 
{V+V-) R (HU. The structure of this first 
order quasiclassical result is identical to that of eq. 
HA6|I . in the latter one should just use the quasi- 
classical form (|29|) for the matrix elements of the 
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E. 



operators ui^2 and replace the coordinates in the 
arguments of the fields by the classical paths, 
V^{Tj,Xj) — > V^{Tj,xin^2m{Tj)). Notc that this 
substitution should be performed neither for the 
field V~ nor for the electron density matrix po 
because no quasiclassical approximation was em- 
ployed with this matrix. Further details are pre- 
sented in Appendix. 



Pauli principle and dephasing in the ground 
state 



It is sometimes argued that electron decoherence at 
T = is impossible because of the Pauli principle: Elec- 
tron at the Fermi surface can neither lose nor gain en- 
ergy, hence, it cannot decohere. Cancellation of "coth" 
and "tanh" terms in the first order perturbation theory 
is considered by some authors as a formal consequence 
of this energy constraint and, on the contrary, indepen- 
dence of on "tanh" terms is interpreted as a sign of 
physical inconsistency of the calculation ( "Pauli principle 
is lost by approximations"). 

Our analysis - which fully accounts for the Pauli prin- 
ciple - does not support the above point of view. Our 
final result, eqs. H4U|) - (|42I) . does depend on the Fermi 
function, however, this dependence enters only into the 
pre-exponent via the matrix elements A"™ (I42II which in 
turn depend on the electron density matrix po. Thus, 
the Pauli principle does not have any significant impact 
on the dephasing process. As we have already explained 
in Ref. |^] and elsewhere, electron dephasing at low T is 
only caused by fluctuations of the field V^. Such fluc- 
tuations are described by the 5i-term in the effective ac- 
tion which is not sensitive to pv at all. In the presence 
of interactions the electron energy fluctuates and it re- 
mains conserved only on average. At the same time elec- 
trons cannot, of course, infinitely decrease their energies. 
Within our formalism such process is prevented by the 
dissipative terms which explicitly depend on py. For in- 
stance, eq. (99) of Ref. [3| demonstrates that electrons 
above the Fermi level decrease energies, however, for en- 
ergies below /i effective "damping" produced by the elec- 
tron bath becomes negative, and the holes are pushed up 
to the Fermi surface. Such processes give rise to the time 
dependence of the pre-exponent contained in the matrix 
elements A'l^ (1^ . 

In the arguments against quantum dephasing at T = 
the Pauli principle is used merely as an energy constraint. 
Therefore such arguments are not specific to Fermi sys- 
tems |22| and can be tested for any quantum particle in- 
teracting with a dissipative quantum environment. It is 
only important to ensure that the whole interacting sys- 
tem "particle-|-environment" is in its true ground state 
at T = 0. One possible way to conduct such a test is to 
study the equilibrium effect of persistent currents (PC) 
for a particle on a ring in the presence of interactions. 
Since nonvanishing PC can only exist in the presence of 



quantum coherence, (partial) suppression of its ampli- 
tude by interactions may signal quantum dephasing. 

Such a problem has recently been investigated by var- 
ious authors and suppression of PC by (lon g: ra i ige) in- 
teractions was demonstrated even at T = [TR llsL l2^ 
(see also ITll'). In p articular, for the model of a diffusive 
electron gas [l3,ll3 one finds that PC gets suppressed by 
interactions exactly in the ground state provided the ring 
perimeter exceeds a finite dephasing length . This 
length turns out to be fully consistent with one found 
from our WL analysis. 

Without going into further details let us briefly ad- 
dress only one point directly related to our discussion. A 
non-perturbative instanton analysis of the problem 
demonstrates that suppression of PC by interactions is 
controlled by the parameter 



Xr. 



(44) 



k=l 



Here A — 3/{8kpP) ^ 1 is the dimensionless interac- 
tion strength, Uk are the Fourier coefficients of the in- 
teraction kernel (a^ ~ (2/7rr) ln(r/fc) for 1 < ^ r and 
flfe « otherwise) and r = R/l ^ 1 with R being the 
ring radius. Provided the parameter (|44|l is large, PC is 
strongly suppressed even at T = 0. The dephasing length 

is derived from the condition Ar ~ 1 which yields |l8l | 

- l/X. 

Up to a numerical prefactor the parameter l|44|) is just 
the instanton action describing tunneling between two 
different topological sectors of the problem. The result 
(I44|l cannot be correctly reproduced within the pertur- 
bation theory. Indeed, let us expand the flux-depending 
part of the free energy and PC to the ffist order in A. 
Then at T = PC is found to be proportional to the 
following combination [T^ 



2 ^ 

fe=i 



kak In 



k + 202 

k — 2(j)a 



(45) 



where —0.5 < (t>x l£ 0.5 is the external flux (normalized 
to the flux quantum) piercing the ring. For small <C 1 
one can expand the logarithm and reduce (|45fl to 



fe=i 



kak 



(46) 



The factors k in the numerator and denominator cancel 
and one is left only with a small correction 2A X]fe=i '^k ^ 
A <C 1. The above cancellation in eq. H46(l obtained 
within the Matsubara technique is to much extent anal- 
ogous to "coth-tanh" cancellation in the real time ap- 
proach. In both cases this cancellation is not com- 
plete, but the remaining term is small and does not give 
the correct answer which can only be obtained by non- 
perturbative means. 

The above example provides yet one more illustration 
of insufficiency of AAV's perturbative approach. For in- 
stance, following AAV's logics one could qualify H44|) as 



9 



"an incorrect perturbative rather than a nonperturba- 
tive" result only because this parameter is proportional 
to the first power of A but does not agree with one derived 
from the perturbation theory H45|) . Proceeding further 
along these lines, one could also "hig hlight" a "mistake" 
in our non-perturbative analysis |l8| . In order to do so, 
one could observe that the same combination A X]fc=i ^^fc 
enters into both non-perturbative (|44|) and perturbative 
(|^ expressions, however the latter also contains the log- 
arithm which is missing in the former. Following the log- 
ics of Ref. one would then be led to conclude that the 
logarithm "is omitted in all orders of perturbation the- 
ory" and the result "is equivalent to calculating only 
a single contribution" A^J^^j^ fcofc. In this way AAV ar- 
rived at their conclusion about missing diagrams 
in our calculation. 

Fortunately, a detailed Monte Carlo simulation pro- 
vides a complete numerical solution for the problem 18]. 
It unambiguously rules out the perturbative result H45|) 
and demonstrates that PC is indeed strongly suppressed 
for Ar ^ 1 even exactly at T = 0, see figs. 1 and 2 of Ref. 
[Tsj l. Similarly, our present results, eqs. lP?Hl - (|I^ . allow 
to discard perturbative calculations of the WL correction 
to conductivity at low temperatures. 



where and are the eigenvalues and the eigen- 

functions of the Hamiltonian Hq. Eq. H48(l was used in 
eq. (43) of and further while constructing the effective 
action. Eq. (|49|l was used in Section 4 and Appendix A of 
Ref. (cf. eq. (54) of that paper) while performing the 
first order perturbative calculation of the conductance. 
With the aid of the form (|49|l in Ref. 5] we have proven 
(partial) cancellation of "coth" and "tanh" terms in the 
first order at T ^ and reproduced the results |3] ■ Also 
the combination l — 2n{p,r) in eqs. (52), (54) and (68) of 
_3] has nothing to do with the "Wigner transform" of the 
density matrix, but is simply equal to tanh(iJo(p, 
This form yields purely real for all paths and Su, = 
for any pair of time-reversed classical paths. 

Let us compare our eqs. (|T7 |) - (B^ with eq. (5) of 
Ref. The latter equation JJ defines an object, pi'4, 
which is neither an operator (cf. our eqs. (|47I48I) ') nor 
the electron density matrix po{r,r')) = (ri|/5o|r2) (cf. 
our eq. H49|l ). We conclude that eq. (5) of Ref. 7] has 
nothing to do with our analysis. Since AAV's claim of 
our "major mistake" and their subsequent critique are 
based on their eq. (5), both this claim and critique can 
be proven irrelevant already by a direct comparison with 
what was actually done in our paper 



VI. REMARKS ON AAV'S CRITIQUE 

The analysis of the previous section not only rules 
out the AAV's claim about vanishing dephasing rate at 
T ^ but also demonstrates that their critique of our 
calculation 0,0] is irrelevant. Nevertheless, for the sake 
of completeness we will reply to both critical points (i) 
and {ii) of Ref. 0. 



A. Density matrix 

In Ref. (31 AAV stated that in eqs. (43) of our paper 
we "replace the density matrix by its "Wigner trans- 
form"", eq. (5) of Ref. '?|. This AAV's statement is not 
correct. The only replacement performed in eqs. (43) 
of m as compared to the exact eqs. (40) of that paper 
(or eqs. H17|l of this paper) is defined by our present eq. 
where 

/5o(p,r) =n(Fo(p,T-)) (47) 

and n(^) = l/[exp(^/r) -I- 1] is the Fermi function. In 
other words, in Ref. 3] we used the following expressions 

1 - 2po(p, r) = tanh ( — ^ I (48) 

and 

(nil - 2po|r2) = ^tanh|^V'.(ri)C(r-2), (49) 



B. Effective action and commutation relations 

In Ref. Ill AAV pointed out that while constructing 
our effective action we disregarded the Poisson brackets 
or, which is the same, the commutation relations between 
the operators po and V~ entering the Hamiltonians H17|l . 
AAV furthermore argued that if one takes care about 
ordering of these operators, one arrives at the effective 
action different from ours. Although the form of this ac- 
tion was not specified, it was claimed in Ref. 7] that the 
term S'r, is not anymore real, but contains an imaginary 
part. According to AAV this imaginary part provides 
nonzero contribution to S'r evaluated on pairs of time 
reversed paths and "in perturbation theory ensures that 
the ultraviolet divergence in iSn cancels that of Si" . 

The latter statement of AAV is false. The correct one 
is just the opposite: It is the reaZ part of S'r that gives rise 
to "tanh" -terms which compensate "coth" -contributions 
in the first order perturbation theory at T = 0. AAV 
seem not to appreciate the fact that the term iSr in 
the exponent of the influence functional and the matrix 
elements generated by this term in the perturbation the- 
ory are different mathematical objects. The perturbative 
contribution from a purely imaginary term zSr can and 
does cancel the contribution from a purely real term Si in 
the first order at T ^ within the Golden rule approx- 
imation. This is a general property not specific to any 
particular calculation. For more information we refer to 
the textbook where the derivation of the perturba- 
tion theory from the influence functional was analyzed in 
details, see eqs. (12-104) to (12-108) of that book. 

As to the commutation relations, everything is, of 
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course, correct with thera in our path integral analysis. 
In order to demonstrate that one should only keep track 
of correct ordering for the operators in the perturbation 
expansion. One way could be to proceed directly with the 
Hamiltonians H17() where ordering is defined uniquely and 
no ambiguity can occur. Alternatively, the full perturba- 
tion theory can be recovered by expanding the influence 
functional in powers of iS^. + Si. In this case one should 
(a) replace the momentum and coordinate variables by 
the operators p ^ p, r ^ f and (b) specify the proper 
way of ordering (fixed by eqs. I17|l 1 in addition to the 
expression for the effective action |25| . 

Furthermore, as it was demonstrated above, in order 
to find Tip at T ~> it is sufficient to correctly derive the 
classical (?i-independent) part of the action only. Obvi- 
ously, while deriving S^'^'^^ there is no need to take care of 
the commutation relations at all. This action is always 
real and is obtained from the quantum Hamiltonian by 
replacing the operators by the corresponding c-number 
functions. For instance, from (|17|) (after the replacement 
Pv{t) Pq) we obtain 



c,(cl) 
•^1,2 



dT 



pi - Uin-,p{x) + eV+(T, x) 

2m 



±-[l-2n{Ho{p,x))]eV-{T,x) 



(50) 



These actions are real and insensitive to the ordering of 
eV~ and 1 — 2po. Hence, the action iSr ^^I obtained from 

s['^2 by averaging over V^, is real as well. On top of 
that, S'r vanishes on pairs of time-reversed paths. Hence, 
it can only contribute to the pre-exponent. The latter 
represents the quantum correction which is sensitive to 
the ordering of the operators. However, this correction 
is formally smaller in the parameter h, and, as we have 
already discussed in Sec. 4, it can never cancel S'^'^^^ as 
long as the latter exceeds Ti. 

In our problem it is not convenient to apply the van 
Vleck formula H29|) directly to the Hamiltonians H12 
(ITTjl . This is because the latter contain the sharp func- 
tion of the electron momentum 1 — 2n{HQ{p,x)) which 
effectively turns fluctuations around the classical paths 
non-Gaussian. This - purely technical - complication is 
circumvented by eqs. (|24|l and the subsequent analysis 
of Sec. 5. This analysis demonstrates that non-Gaussian 
fluctuations give rise to the pre-exponential factors Ai^2 
in the expression for the conductivity (|40|1 . We have 
proven in Sec. 5 that these factors are irrelevant for de- 
phasing because they do not grow at long times and, 
hence, cannot cancel the term Si. 



VII. SUMMARY 

In summary, we have derived a complete expression for 
the weak localization correction to the conductivity of 
a disordered conductor in the presence electron-electron 
interactions. Our analysis has been carried out within 



the standard model for an interacting electron gas in dis- 
ordered conductors with no approximations beyond the 
accuracy of the definition of the WL correction. In partic- 
ular, interactions have been treated non-perturbatively, 
no quasiclassical approximation for the electron density 
matrix has been employed and no disorder averaging has 
been performed at all. We have fully confirmed our ear- 
lier results H, H, IE IB] and extended them by explic- 
itly taking into account quantum fluctuations around the 
classical paths for interacting electrons. We have proven 
that such fluctuations, while practically irrelevant for the 
calculation of Tip, do contribute to the Cooperon dynam- 
ics at short times causing, for instance, partial cancella- 
tion of the well known "coth" and "tanh" terms in the 
first order perturbation theory. We have also demon- 
strated the failure of a perturbative calculation 7] in 
the problem of quantum dephasing of electrons at low 
temperatures. Finally we have refuted AAV's critique of 
our previous calculation 0,01 observing that (i) Poisson 
brackets are irrelevant for the problem of electron dephas- 
ing by interactions and {ii) no "Wigner transform" of the 
electron density matrix was performed in our derivation. 
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Note added. After this paper had alreadybeen submit- 
ted there appeared an independent work |2g addressing 
the same issue. In this work von Delft (vD) has success- 
fully re-derived our influence functional for interacting 
electrons Q and argued that (a) our approach "properly 
incorporates the Pauli principle" and (b) "the standard 
Keldysh diagrammatic expressions for the self energy of 
the Cooperon can be obtained from iSn -\- 6*1" , i.e. from 
our influence functional. Thus, it is now verified not 
only in our work but also independently by other authors 
|26l l27j that our path integral result H10|l Q contains all 
RPA diagrams to all orders in the electron-electron in- 
teraction. 

The observations (a) and (b) are important because 
they allow to restrict the whole discussion to just one - 
purely mathematical - issue, i.e. how to correctly evalu- 
ate the path integral 1)10(1 . We believe that the analysis 
presented in Sec. 5 of this paper should eliminate all 
doubts ^23| concerning the role of S'^-terms for quantum 
dephasing of electrons at low temperatures. This analy- 
sis, for instance, involves none of the approximations de- 
noted in as (i), (ii) and (iii). In particular, it rules out 
vD's conjecture that within our approach we "neglect all 
the diagrams of Fig. 2b" of 26] . Quite on the contrary, 
our final result, eqs. H40|l - H42l) . explicitly accounts for all 
these diagrams (giving rise to the "tanh" -contribution 
(|A9|) in the first order) as well as for infinitely many di- 
agrams of all higher orders not presented in Ref. [2g. 
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Several additional comments are in order: (i) The 
statement |^ that the first order perturbative result 
contains no ultraviolet (UV) divergences is explicitly in- 
correct for Id and 2d systems, see, e.g., eq. (70) of Ref. 
0. (m) We disagree with vD's conjecture (see the foot- 
note 16 in Ref. |2^) that diagrams with crossed and 
overlapping interaction lines can be neglected |28^. (iii) 
Further evidence that the above conjecture is problem- 
atic is provided by the results which we also address 
in Sec. 5E In that problem the first order diagrams 
yield negligible contribution (|46|l and, hence, the correct 
result (|44|l is dominated by all the remaining diagrams 
with crossed and overlapping interaction lines, (iv) The 
argument presented by eqs. (10-11) of 26] is inconclu- 
sive, since it is based on an improper application of qua- 
siclassical methods to Fermi systems 30] . If one expands 
the electron density matrix po in powers of h one indeed 
gets a series of terms diverging at T ^ 0. However, this 
observation can only imply that the Taylor expansion of 
the step function (i.e. the Fermi function at T ^ and 
energies close to n) is essentially useless. A much more 
useful strategy is to retain po (which is, of course, always 
finite and not large) in its full quantum mechanical form 
and to apply quasiclassics only to those matrix elements 
which do not contain pQ. This strategy was implemented 
in Sec. 5. {v) In contrast to vD's conjecture in the foot- 
note 23 of 26] our result (^UJ-lESl does not diverge at 



APPENDIX A 

Here we will summarize several expressions used in our 
calculation and present some perturbative in the interac- 
tion results for the conductivity. 

Averaging over the fluctuating fields implies calcu- 
lating the double path integral 



{...)v^ J VV+VV-{...)exp^-SEM[V^]j ■ (Al) 

Here Sem can be understood as a formally exact effective 
action for the Hubbard-Stratonovich fields , see, e.g., 
eq. (10) of Ref. 3]. For the situation discussed here it is 
sufficient to expand Sem to the second order in V. Then 
one finds 



2 J (27r)' 



-V-{-K)- 



coth — V-{K)Ak2) 
47r 2T ^ ''^ ' 



where K — {uj,k). This action allows to determine the 
correlation functions 



{V+{t,r)V+{0,Q))^M{t,r) 



-in 



= ^ I T-^ Im ( ) coth — e^'-^^, (A3) 

(27r)4 \k^e{K) ' "^^ ' 



2T 



{V+{t,r)V-{0,0)) =inR{t,r) 



— ih 



{V- {t,r)V- (0,0)) ^0. 



d'^K At: 



iKX 



(27r)4 k^eiK) 



,(A4) 
(A5) 



Here we have defined X — {t, r), and KX — uot — kr. Em- 
ploying the action 1A2|I is equivalent to describing the 
electron-electron interaction within the random phase 
approximation. 

Let us expand the results H14I24() to the first non- 
vanishing order in the interaction (second order in V^) 
and replace pv{t') = Po- Then we obtain 



'3h 



dt 



' (tr {j{x)uQ{t,t')[x,pQ\uQ{t' ,t) 



dt' / dri I dr2lm 



tr[j{x)uo{t,Ti) 



2e^ 
3^ 



xy+(ri)uo(Ti,r2) <; V+{t2) + ^{l ~ 2po)y-(T2) 

XU(3{T2,t')[x,PQ]uQ{t\t) 



+iAj{x)uo{t,T2) V+{t2) + -(1 - 2po)V-{T2) 



Uo{T2,t')[x,PQ]uo{t' ,Tl)V^(Tl)U(^{Tl,t) 



(A6) 



After averaging over the latter expression coincides 
with the result derived in Ref. 3 . The term proportional 
to coth|y emerges from the average (y^ [ti)V'^ {T2)) ^ 
while tanh ^2t" a-ppears from the combination 1 — 2pQ. 
Further details can be found in Ref. 0- 

Note, that eq. HA6|I represents the exact first order 
result obtained without any evaluation of the path inte- 
grals. Eq. (|40|l is valid to all orders in the interaction, 
but it was derived by evaluating the path integrals in the 
quasiclassical limit kpl ^ 1. Let us expand J"™ H41|) to 
the first order in 5*1 and l|42|l to the first order in 

Vq. Then with the aid of eq. 1401) we reproduce eq. (|A6() 
with trivial modifications as described towards the end 
of Sec. 5D. 

Let us now identically transform our first order quasi- 
classical results for <5(Tj^^ to a somewhat different form. 
For that purpose let us express the electron density ma- 
trix Po as follows: 



{zi\Pq\^2) = Po(zi, 2:2) 

+ CXJ 



ds2 



iT 



Uo{-S2,Zi,Z2),iA7) 



where Uq{—s,zi,Z2) = (2:i|wo(0, s)|z2) is the matrix ele- 
ment of the evolution operator for non-interacting elec- 
trons. In addition we note that the functions I{t, r) (|A3|) 
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and R{t, r) ljA4p are related to each other by means of 
the identity 

r T coth 

I(X) = J dsi —J^[R(t^Si,r)+R{-t-s,,r)]. 

I 



Eqs. HA8|I and (jA9|l have a very similar structure. These 
two expressions are, however, not fully identical even at 
T —f and, hence, they do not cancel exactly in the first 
order result H43I) . see also Sec. 4 of Ref. 

Further evaluation of eqs. (|A8() and HA9|I makes little 
sense because the first order perturbation theory can- 
not provide any useful information about the electron 
dephasing time at low temperatures. Nevertheless, the 
above expressions are of a certain interest, since they help 
to illustrate the relation between perturbative and non- 
perturbative results at the stage when the quasiclassical 
approximation has already been performed. We observe, 
for instance, that all the first order terms, both "coth" 
and "tanh" contributions, are fully reproduced from our 
path integral analysis. Another observation concerns the 



Making use of the above identities we arrive at the fol- 
lowing expressions for (5(t|\^: 



(A8) 



(A9) 



relation between the quasiclassical paths emerging from 
the path integrals and those entering the first order re- 
sults for SalP . The WL correction to the conductivity 
is defined on pairs of time-reversed path, and only such 
paths (plus fluctuations around them) are relevant for 
the path integral analysis of this quantity. Of course, the 
same paths enter if the general result is expanded to the 
first order in the interaction before the transformation 
(|A7() . However, after this transformation there appear 
additional matrix elements uq of the electron evolution 
operator. Proceeding quasiclassically, one can evaluate 
these matrix elements by means of the van Vleck formula 
(|25J), i.e. to write 

uo(i'-T2-si,?/i,r) cxe^^o"(*'-^2-si,o,2,i.r) ^^^q^ 



and 
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X < I dri J dT2 
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X (2/1 - y2)uo{~'S2,yi,y2) 



X (1) 



/ dt' / dsids2 / dyidy2dzdr 



12mh 
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dxidyi 



dx2dy2 
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I 
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and similarly for other matrix elements. Substituting uq 
in the form HA10|1 into eq. HA9|I one can interpret the 
result in terms of the electron motion along additional 
classical paths ifc(s), say, first from yi to r and then 
from r to z (some of these paths violate the requirement 
of causality, see Fig. 3 of Ref. and related discussion) . 



This could in turn create an illusion that these additional 
paths are missing in the path integral formulation. The 
above analysis clearly indicates the origin of such an illu- 
sion. It also demonstrates that - in contrast to Ref. Q - 
the whole issue has nothing to do with disorder averaging 
which is not performed here at all. 
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